Photon heat transport in low-dimensional nanostructures 
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At low temperatures when the phonon modes are effectively frozen, photon transport is the 
dominating mechanism of thermal relaxation in metallic systems. Starting from a microscopic 
many-body Hamiltonian, we develop a nonequilibrium Green's function method to study energy 
transport by photons in nanostructures. A formally exact expression for the energy current between 
a metallic island and a one-dimensional electromagnetic field is obtained. From this expression we 
derive the quantized thermal conductance as well as show how the results can be generalized to 
nonequilibrium situations. Generally, the frequency- dependent current noise of the island electrons 
determines the energy transfer rate. 



General physical and information-theoretic arguments 
imply that there is a fundamental limit Gq — ir 2 k^T/3h 
to the thermal conductance of a single channel^, in- 
dependent of the nature of the conduction mechanism. 
Particularly, Gq should be independent of the disper- 
sion relation and quantum statistics of carrier particles. 2 
The few-channel heat conductance is particularly rele- 
vant in low-dimensional nanostructures where the chan- 
nel number is naturally low. The quantized thermal con- 
ductance has been experimentally verified for electrons, 
phonons and recently in photon transport between metal- 
lic islands!^£ If electron transport is restricted and the 
system is at very low temperature so that the phonon 
modes appear frozen, the dominant thermal relaxation 
process is photon transport 

In this Letter we study energy transport by photons 
between a metallic island and a one-dimensional electro- 
magnetic field supported by a transmission line. The 
latter mimics the effect of the external leads and con- 
nectors on the island. Our aim is to give the photon 
transport a microscopic basis as well as to study nonequi- 
librium processes. Applying Green's function methods to 
a microscopic model, we obtain a formally exact expres- 
sion for the energy current. We show that frequency- 
dependent current noise determines the characteristics 
of the transport process, thus providing a close connec- 
tion between the electrons and the photons. We derive 
an expression for the heat flow between the field and 
the island and verify that the maximum value of thermal 
conductance in the system is Gq . This provides a micro- 
scopic description of the recent experiment on electron- 
photon coupling^ We consider also a many-channel case 
where the electron system is connected to several trans- 
mission lines. The energy current formula allows us to 
study situations where the island is driven to a nonequi- 
librium state and examine how electron shot noise al- 
ters the energy transport. We show that, due to shot 
noise, part of Joule heat flows to the photons. Our results 
have relevance in determining the electron temperature of 
driven mesoscopic systems^, but they are also important 
in photon-based solid-state applications, such as cavity 
QED and its quantum information realizations^^ 

The studied model consists of a small metallic island 
coupled to a parallel strip transmission line, see Fig. [T] 



FIG. 1: (Color online) Metallic island (blue) coupled to the 
electromagnetic field of a transmission line (lines around the 
yellow region). The voltage between the strips is Vtl- 



The transmission line acts as a waveguide supporting 
a one-dimensional electromagnetic field. In contrast to 
three-dimensional waveguides, the parallel strip line field 
has only one allowed field polarization. Therefore it cor- 
responds to a single transport channel. There is no direct 
electrical connection between the electrons on the island 
and those in the transmission line strips, only the field 
couples to the electrons. The total Hamiltonian of the 



system is H = H e + H 1 



iJ e _ 7 , where 
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In the following, we do not have to specify the terms in 
the electron Hamiltonian H e in more detail. The trans- 
mission line is characterized by its length L, distance be- 
tween the parallel strips L z , and inductance I and capac- 
itance c per unit length. Operator Vtl = J2j Tj{aj + fit) 
is the voltage operator at the end of the line. The inte- 
gration in -f/c-7 is restricted to the region between the 
parallel strips in the case the electron system extends 
beyond that. The island is assumed to be much smaller 
than the photon wavelength at the relevant frequencies, 
so the position dependence of the voltage operator in 
the interaction term can be neglected. The field opera- 
tors $ (r) , (r) and creation and annihilation operators 
aj , fit satisfy canonical fermion and boson commutation 
relations. Constants Wj = jirv/L (j is a positive integer), 
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Tj = y/ftiOj/Lc and g = e/L z can be found by quantizing 
the line field.— The wave velocity v in the transmission 
line is given by v = 1/ Vic. 

The electron Hamiltonian H c does not commute with 
the total Hamiltonian H, thus there is an energy flow 
between the island and the field. This energy flow into 
the island is characterized by a current Jq defined as 



J Q = (H e ) = -2- l^j &{r)pA(r)drV TL 



(4) 



The notation (•) stands for averaging over a density ma- 
trix of the total system. We calculate averages over 
noncquilibrium states where subsystems have a temper- 
ature gradient or electron system is subjected to a finite 
voltage. To simplify expressions, the following short- 
hand notations are introduced: P = J (r)p z ^>(r)dr, 
Z = J &(r)zi>(r)dr and / = -^P- The quantities are 

related by Z = P/m. The current (|4]) can be written as 



JqV) 



2fJ 



(5) 



The energy transport problem is reduced to finding the 
"lesser" Green's function Gf{t,t'). We will concentrate 
on a steady-state situation where Gf(t,t') — Gf{t — t'). 

The contour-ordered Green's function Gj(r, r') can be 
derived by the equation-of-motion technique Let us 
first consider the time-ordered Green's function G^it—t') 
at T — 0. By differentiation and applying Heisenberg's 
equation of motion we obtain 



{id v -u> j )G){t-t') = y -T j {P{t)Z^)Y 



(6) 



The expression in brackets on the left hand side of Eq. © 
can be interpreted as an inverse Green's function ■~ 1 
of a free photon field. Thus Eq. © can be solved by 
integration, yielding 

G'At-t') = 9 -T 3 [ dii(P(t)Z(ti))*Dj(ti -0- 



Using analytical continuation rules known as the Lan- 
greth's theorem^ we obtain 

Gf(t - = 9 -T 3 J dtt [(PW&WDffa - t')+ 

+(P(t)Z(t 1 ))<D?(t 1 -t') 

Superscripts r , a and < stand for "retarded", "advanced" 
and " lesser" . For later purposes it is convenient to define 
the Fourier transform 

G<M = lgT s [(PZy(u)D<{u) + (P2)<(u)D5(u) 
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(ii) r (u)Df(u) + -(ii)<(u>)D?(u) 
j .i j 



Now Eq. ((5l) yields for the steady-state current 

did 
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For a transmission line much longer than Kitv/QcbT), the 
sum over the field modes can be replaced by integration 

according to J2j = k fo° ^ = mi Ia° ^ ne curren t 
takes the form 



Jn = Ke / dLOjfujjj I —x 

n J J 2n 



two 



(8) 



where Zq = sfljc is the characteristic impedance of the 
transmission line. 

The photon Green's functions at a finite tempera- 



ture can be written as D"> = — 27rin 7 (u)5(aj — u>j) and 
m = 1 — 

J u—ujj—ir} 



iri5(ui — uij) + P u ^ , where n 7 (u) is 
the Bose distribution. Inserting these expressions into 
Eq. |8]| gives 



Jq = 2Zq 



du> 
2^ 



2Re(ii) r (tu)n 7 (uj) - (ii)<{Lu) . (9) 



The correlators on the right hand side of Eq. ([9]) 
can be expressed in terms of the noise power Si — 



I. 
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(I(t)I(0))dt as (II) < (lo) = Si(-w) and 



Re(II) r (uj) = - S/(-o;)). Expression ® is an 

exact formula for the energy flow and valid even when 
the electron system is out of equilibrium. However, it 
contains the exact current-current correlation functions 
of the metallic island in the presence of the field. In 
equilibrium, these are related to conductance through 
the Fluctuation-Dissipation theorem and the Kubo for- 
mula. From formal point of view, the exact expression for 
noise power determines the energy exchange process com- 
pletely. In the weak-coupling limit (the lowest order in 
electron-photon coupling), one can neglect the field and 
use the bare island correlators. On physical grounds one 
expects that the maximum energy transport is achieved 
when the coupling is strong, thus a more accurate treat- 
ment of the electron-photon interaction is desirable. 

The quadratic form of H 7 and the linear coupling term 
H c -ph, together with the density of states of a long trans- 
mission line is precisely a Caldeira-Leggett representation 
of an ohmic losspii Solution to the equation of motion of 
the voltage operator is Vrh(t) — VShi^) + Zal{t), where 
Vj! L (t) is the solution in the absence of the island and 
I(t) is the current flowing in the electron system. This 
notion microscopically motivates the circuit approxima- 
tion, where the transmission line can be thought as a 
resistor in series with the metallic island, see Fig. [5] a). 



In the circuit description correlation functions can be 
calculated by the Langevin approach, which allows us 
to relate the current fluctuations in the presence of the 
environment to bare quantities^ This yields 



\l + G e (u)Z \ 



(10) 



where (II) e (u>) and G e (ui) are the current-current cor- 
relation function and conductance of the island in the 
absence of the electromagnetic field. According to the 
Kubo formula for conductance, the real part of the re- 
tarded function appearing in Eq. (J9]) is related to the 
conductance as Re[G e (w)] = Re[(II) r e (uj)}/{hiu). With 
approximation (|10|) . we then get 



Jq 



duj 



2Zn 



2tt |1 
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G e (cu)Z \ 2 
Sf(-u;)]n» 



where Sf(ui) is the noise power for the isolated electron 
system. 

In (quasi)equilibrium the correlation functions are re- 
lated through a variant of Fluctuation-Dissipation for- 
mula Sj(lo) = 2Re[G e (uj)]hujn e (uj), where n c {oj) is the 
Bose distribution function at the electron temperature. 
Thus for this case 



Jq 



4Z R C 



(Re + z ) 2 J 



dujhuj 
2tt 



n 7 (w) - n e (u>)] , (12) 



when the island is assumed resistive R e = 1/G e . Result 
(fT2"| agrees with the one stated in Ref. @ for heat flow 
between two resistors. After integration Eq. (p~2|) gives 



Q 



r^(T 2 -T 2 ), 



4Z R C 



6h y ^ ~ c " ' (R C + Z Q ) 2 
At small temperature difference this is just 
J Q =rG Q AT, 



(13) 



(14) 



where Gq = ir 2 kgT/3h is the universal quantum of heat 
conductance^. Thus, when R e = Zq and thus r = 1, the 
maximum one-channel heat transfer is achieved. In the 
weak-coupling limit, where the exact correlation func- 
tions in are replaced by bare correlators (II)o(oj), 
one recovers (fl4)) with the prefactor Z^R C /(R C + Z ) 2 
replaced by Z /R e - Physically the weak coupling result 
follows from the impedance mismatch Zq <C R e - 

The above discussion can be generalized to incorporate 
several photon channels realized by coupling the elec- 
tron system to, say, N transmission lines, as in Fig. [5] 
Suppose that each transmission line is described by a 
Hamiltonian of the form ^ with the coupling (J3J) corre- 
sponding to the situation in Fig. [2]b). The theoretical 
maximum heat conductance for N independent channels 
is N x Gq, but it is not achieved in this case. An added 
transmission line does not simply add an independent 
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FIG. 2: (Color online) a) Circuit picture corresponding to 
the studied system. Transmission line acts as a series resistor 
to the island. In b) the parallel transmission lines couple to 
vertical current noise and in c) the lines couple to vertical and 
horizontal noise. The maximum heat conductance in b) is Gq 
and in c) is 2Gq. The number of parallel transmission lines 
is irrelevant for the maximum heat conductance. In d) the 
island contains a short conductor and is externally biased. 



photon channel because it also effectively acts as a se- 
ries resistor in the coupling direction. Thus it suppresses 
current fluctuations and affects the emitted energy in all 
channels. The heat flow (fl4|) for multiple channels is 



J Q - /J G C 



mZi 



(15) 



where Zi is the characteristic impedance, AXi the tem- 
perature difference and Ri the electron resistance asso- 
ciated with line i. When all the transmission line fields 
are at the same temperature, the maximum heat conduc- 
tance given by Eq. (fl~5|) is still Gq. Thus, adding parallel 
lines does not increase this maximum. However, cou- 
pling the island to perpendicular transmission lines, as 
in Fig.[2]c), opens up an independent transport channel. 
The difference in b) and c) is that the lines in perpendic- 
ular directions couple to different current components. 
The flow Jq is then a sum of two terms of the form (|T5|) 
and the maximum heat conductance is 2Gq. Similarly, 
coupling to the remaining orthogonal direction yields the 
maximum heat conductance 3Gq. 

Next we consider a case where the island contains a 
short contact which is externally biased by potential dif- 
ference, see Fig. [2]d). Supposing that the electron trans- 
port is coherent and neglecting interaction effects, cur- 
rent noise contains both the equilibrium and shot noise 
and can be written 



S I (co)=G o yT m (l-T m ) 



m 

-eV + Huj 

I _ e Ji(-hui+eV) 



eV + hw 



Gc 



e 0{-huj-eV) 

r^ T o 2fuv 



(16) 



/ . m -i 
— 1 — 



where Go = e 2 /h, V is the bias voltage and T m is the 
transmission eigenvalue of channel m. The sums of trans- 
mission eigenvalues extend over the channel index and 
the spin. Inserting expression (| 16[) to the general for- 
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mula ([9]) we discover 



p 



G C (T^-T^--F 2 GV' 



(17) 



where G = G J2 rn T m — X/R e is the island conduc- 
tance, F 2 = J2m T m( 1 - T m)/Y, m T m the Fano factor, 
and Cq = 7r 2 fcg/3e 2 the Lorenz number. The last term in 
Jq corresponds to the increased emission by shot noise. 
The frequency dependence in Eq. (fT6]) is solely due to 
the Fermi distribution and the emitted energy due to 
shot noise shows only dependence on the bias voltage. 
Expression (|16[) is valid only for low frequencies; gener- 
ally Si(u) probes the intrinsic (inverse) time scales of 
the conductor such as the time of flight and the charge 
relaxation time. This is shown in Fig. [3] where we have 
used the noise and conductance of an interacting chaotic 
cavity^ to numerically compute the energy flow. 

Assume that the island is biased using superconduct- 
ing wires with contact conductances much higher than G. 
Such a setup provides thermal insulation of the island 5 
while the voltage still drops across the contact. The final 
temperature T e of the island can be obtained from a heat 
balance equation where the Joule heating Jq from the 
voltage source is balanced by heat flow to photons and 
phonons as J^ + J? h = GV 2 +J^+J:n[T^T^} = 0, 
where S is the electron-phonon coupling constant! and f2 
is the volume of the island. There is a crossover temper- 
ature T cr = (rG £o/(2£^)) 1/3 below which the photon 
transport is the dominant process. For example, with 
the parameters of Ref. H T cr would be roughly 140 mK - 
for smaller objects such as carbon nanotubes, it could be 
made larger at least by one order of magnitude. Much 
below T cr the final electron temperature is 




and above the crossover it is 



T, 



tItI 



(18) 



G U--5- ™ ■ (19) 



F2T\ V 
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In both limits the Joule heating is reduced by the factor 
(1 — F 2 r /2) because a fraction of it flows to photons. In 
the case of an ideally matched (r = 1) tunnel junction 



(F 2 = 1), exactly half of the Joule heat goes to photons. 
In experiments the parameters V , T 7 and r can be 
varied in situ to investigate the photon transport contri- 
bution, as demonstrated in Ref. 5 for r, T p h and T 7 . 

In conclusion, we studied a microscopic model of pho- 
ton transport in nanostructures using a Green's-function 
method and derived a general expression for the energy 
flow between a metallic island and a transmission line 
field. We showed how electron and photon transport are 
related through frequency-dependent current noise. We 
demonstrated efficiency of the energy flow formula by de- 



riving quantised nhoton heat conductance and studying 




FIG. 3: (Color online) Energy flow from a symmetric chaotic 
cavity (Nl = Nr = 1) as a function of voltage. We as- 
sumed that the cavity conductance at zero frequency G e (0) is 
matched to Zq . The different curves correspond to different 
cavity charge relaxation times r, r = (solid), r = 0.1/1/3 
(dashed), t — hp (dotted), r = 10hf3 (dash-dotted). Inset 
shows the r/ft/3-dependence of the heat flow from the cav- 
ity (V — 0), the horizontal dashed line corresponds to the 
heat flow from an ideally matched ohmic resistor without fre- 
quency dependence. When t/H/3 <C 1 the heat flow is close to 
the theoretical maximum and settles to a lower value as the 
fraction increases. 



effects of electron shot noise to photon transport. We 
propose to measure the shot-noise effect illustrated as 
the voltage-dependent term in Eq. (|17|) by modifying the 
setup in Ref. lfj| to include a small mesoscopic junction. 
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